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Abstract 



A polygon representing transportation network is given, together with a point p 
£N| in its interior. We aim to extend the network by inserting a line segment, called a 

feed-link, which connects p to the boundary of the polygon. Geometric dilation of 
some point q on the boundary is the ratio between the length of the shortest path 
from p to q through the extended network and their Euclidean distance. The utility of 
a feed-link is inversely proportional to the maximal dilation over all boundary points. 
We give a linear time algorithm for computing the feed-link with the minimum overall 
dilation, thus improving upon the previously known algorithm of complexity close to 
£S| 0(n\ogn). 

o , . 

1 Introduction 

CO 

i— i Given a polygon and a point p inside it, we want to connect p to the boundary of the 

polygon using a single line segment. Any such connection is called a feed-link. The optimal 
j> feed-link is the one that minimizes detour factor from point p to any other point on the 

boundary. 

G\ 

The problem of finding the optimal feed-link is solved in 1 1 j in 0{\j(n) logn) time, where 
\j{n) is a slightly superlinear function. We make an improvement by presenting an 0{n) 
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time algorithm. 

Although the problem statement assumes that p lies inside the polygon, all the calculation 
work out exactly the same for an arbitrary point p in the plane, and the same result is 
obtained in that more general setting. 



2 Notation and problem statement 



Polygon with boundary P is defined by vertices Vo, v%, . . . , v n -\. The point p in its interior 
is given. 

Feed-link is any line segment pq, connecting p with some point q G P on the boundary. 
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Given points q,r £ P dilation of r via q is denned as 

\pq\ + dist(g, r) 



S q (r) 



\pr\ 



where dist(g, r) is a distance between q and r by going over the polygon's boundary, and 
\ab\ is the Euclidean distance between points a and b, see Figure [TJ 




Figure 1: The concept of dilation. 



Given point q £ P, dilation via q is defined as 



S„ = max^ g (r). 

The problem of finding the optimal feed-link is to find q such that 8 q is minimized. 
2.1 Left and right dilation 

We split the concept of dilation into the left and right dilation. 

Given two points a,b G P, P[a, b] is the portion of the P obtained by going from a to b 
around the polygon in the positive direction, including points a and b. Let (i(a, b) be the 
length of P[a,b] and fJ.(P) is the length (perimeter) of P. 

For given a £ P, let a' be the point on P, different from a, for which /i(a, a') = n(af , a) = 
fi(P)/2, see Figure [2] By P + [a] we denote P[a,a'], and by -P~[a] we denote P[a',a\. 
Obviously, P+[a) U P~[a] = P and P + [a] n P~[a) = {a, a'}. 



P+\a 




p-[a] 

Figure 2: Left and right portion of P from the point a. 

Given point q £ P and r G P + [q] , Ze/t dilation of r via q is defined as 

\pq\ + r) 



|pr| 



Given point q £ P and r £ P [q], right dilation of r via q is defined as 

= jggj + Kr^g) 
9 |pr| 
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When measuring dist(g, r) we take the shortest path from q to r over P, and that whole 
path must lie entirely either in P + [q] or P~ [q] . This allows us to express the dilation of r 
via q using left and right dilations of r via q 



S q (r) 



5+{r) if reP+[q} 
6-(r) ifr€P-[g] 



Given point q £ P, left dilation via q is defined as 5^ = max rg p+[ (? ] 5£(r), and right 
dilation via q as 5~7 = max rg p-[ 9 ] 57 (r). Finally, the dilation via q can be expressed as 

5 q = max(<5j~, 5~) = max5 9 (r). (1) 

In the next two sections we will be concerned only about left dilations; analogous result will 
hold for right dilation. In Section [5] we will show how to combine our findings about left 
and right dilation to compute the answer to the stated question. Here we stop introducing 
notations with + and — in superscripts denoting left and right dilations, and consider each 
new notation to be related with left dilation, if applicable. 



3 Another view of the problem 



We parametrize points on P by defining P(t), t 6 1, to be the point on P for which 
/j,(va,P(t)) = t (mod /u(P)), see Figure [3j 




Figure 3: Parametrization of P. 



The plot of the function h(t) = \pP(t)\ is an infinite sequence of hyperbola segments joined 
at their endpoints, where (kn + r)-th hyperbola segment corresponds to the r-th side of 
P, for r G {0, 1, . . . , n — 1} and k G Z, see Figure |4j For each i = kn + r, hyperbola hi 

is of the form hi(t) = \J {t — rm) 2 + df, where mkn+r = fn r + kfi(P), and dk n +r = d r , d r 
being the distance between p and the line containing r-th side of P. The left endpoint 
of i-th hyperbola segment is Ei = (ej,/i(ej)), where e« = k/j,(P) + n(v r ), and the right 
endpoint is at (ej + i, /i(ej + i)), so < e$+i. We will consider that each hyperbola segment 
contains its left endpoint, but not the right endpoint. Let H(t) = (t,h(t)). We denote 
i-th hyperbola segment with Hi- The plot is, obviously, periodic, with the period of [J-(P), 
that is, H(t) = H(t + kfi(P)). 

Let o(t) = t - h(t), and 0(t) = (o(t),0). We also define o;(i) = t - hi(t), and Oi(t) = 
(oi(t), 0). For given t q and t r , such that t q < t r < t q + /j,(P)/2, the slope of the line passing 
through points 0(t q ) and H(t r ) is 

a t )= h{tr) = lpP{tr)l = 1 
u ' r) t r -t q + h(t q ) v(p(t q ),p(t r )) + \pp(t q )\ 5+ t) (p(t r )y 
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t 

Figure 4: The plot of h(t). 



hence the slope between 0(t q ) and H(t r ) equals the inverse of the left dilation of P{t r ) 
via P(t q ), see Figure [5j 




r 



Figure 5: Dilation and slope relation. 

The lowest slope s(t q ) from 0{t q ) to H(t r ) among all t r £ [t q , t q + fj,(P)/2] thus equals the 
inverse of the left dilation via P(t q ): 

s(t ) = min s(t g ,t r ) = min — r — - — — 

! ! (?) (2) 

max r . eP+[P(t9)] 5+ M (r) 5+^' 

Obviously, s(i) G (0, 1] because it is strictly positive and s(t) < s(t,t) = 1. 
Lemma 1. For any two distinct values of t± and t2, 

|(M*2)-M*l))/(*2-*l)| < I- 

Proof. Function h{t) is continuous and, as a union of countably many segments of hy- 
perbolas with first derivatives less or equal one in absolute value, is differentiable almost 
everywhere having < 1 for each i £ l\{e, : i £ M}. □ 

So far, s{t q ) was defined as minimum only among slopes s(t q , t r ) where t r G [t q , t q +fj>(P)/2]. 
However, this range can be extended since from Lemma [T] follows that s(t q ,t r ) cannot be 
less than 1 when t r £ [o(t q ),t q ]. Hence, 

s(t a ) = min s(t a ,t r ). 

t r e[o(t q ),t q +^P)/2] 

4 Sliding lever algorithm 

Consider the line segment having slope s(t), with one endpoint at 0(t) and the other at 
(t + jU(P)/2, s(t)(n(P)/2 + h(t))). Let us call that line segment the lever for t (Figure [6}. 
Note that the lever only touches the plot, never intersecting it properly. 
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Let C(t) = H(c(t)) = (c(t) , h(c(t))) be the leftmost point in which the lever for t touches 
the plot. Then, c(t) £ [o(t),t + fi(P)/2] is the lowest value for which s(t) = s(t,c(t)). 
From the original perspective, this means that left dilation via P(t) reaches maximum at 
P(c(t)). 




0{t) t t + n/2 



Figure 6: Lever. 



We continuously decrease parameter t and observe what is happening with the lever for 
t. Decreasing t corresponds to "dragging" the lever strictly in the leftward direction, 
because when t is decreasing, the following monotonicity lemma states that o(t) and c(t) 
are decreasing as well. 

Lemma 2. t\ < t 2 implies o(ti) < o{t 2 ) and c(t%) < cfa) 



Proof. Suppose t\ <t2- From Lemma [TJ 



(h(t 2 ) - M*i))/(*2 - h) < 1 

t x - h{ti) <t 2 - h{t 2 ) 
o{ti) < o(t 2 ). 



To show that c(t{) < c(t 2 ), assume the opposite, that t\ < t 2 and c(ti) > c(t 2 ). Then, 
h < t 2 < c{t 2 ) < c(tx) <h + fi(P)/2, and c{h) G [t 2 ,t 2 + n{P)/2\. 

If line segments 0(ti)C(ti) and 0{t 2 )C(t 2 ) do not intersect, then 0(t 2 )C{t 2 ) lies com- 
pletely under the 0(ti)C(ti) and, since 0(t 2 )C(t 2 ) touches the plot, the plot must inter- 
sect 0{ti)C{ti) in some point left of ti, which is impossible since C{t\) is the leftmost 
point where the lever for t\ touches the plot. 

Otherwise, suppose 0(t\)C[t{) and 0{t 2 )C{t 2 ) intersect. 

s(i 2) c(ti)) < s(t 2 ,c(t 2 )) = min s(t 2 , t) < s(t 2 , c(tt)), 

te[t2,i2+M(^)/2] 

which is a contradiction. □ 



We now introduce some functions to be used later. 

Let Cj(o) = H(cj(o)) = (cj(o), h(cj(o))) be the contact point of hyperbola hj and its 
tangent through the point (o, 0), o < rrij. It can be calculated by solving the equation 
h'j(cj(o)) = hj(cj(o))/ (cj(o) — o), which results in 

dh + m 2 — rrijO 

c» = " -• (3) 

m,- - o 
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For example, if we know that the lever for t is tangent to the hyperbola segment Hj, we 
then know that it is touching it at the point Cj(o(t)). 

Next, let Sj(o) be the slope of that tangent, the line through (o, 0) and Cj{o). 




4.1 States 



We define different lever states depending on where t and c(t) are. 

When t G [e*, ei+i) and c(t) £ [ej, ej+i), we say that the lever for t is in state (£, j). There 
are three possible types of positions defined by the manner in which the lever for t touches 
the plot: 



• (hj) '■ c (t) < t + and the lever is the tangent to %j. The lever is sliding 

to the left along the hj maintaining the tangency, thus continuously decreasing the 
slope. 




• {i,j)y : c(t) = t + fi(P)/2. Point C(t) is the right endpoint of the lever. It is the 
only point where lever touches the plot. The lever is moving to the left while keeping 
its right endpoint on hj. 




• {i,j) v : c(t) < t + n(P)/2 and the lever is passing through the point H(ej), the 
endpoint between hyperbola segments Hj—i and Hj. This situation occurs only 
if rrij-i > rrij. The two neighbouring hyperbola segments then form a "wedge" 
pointing downwards, and the lever is sliding to the left while maintaining the contact 
with the tip of that wedge, thus continuously decreasing the slope. 
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4.2 Events 



During the process of dragging the lever towards the left, the lever state changes at certain 
moments. We call such events state transition events. If the current value of t, denoted by 
t c , and the current state are known, the following event can be determined by maintaining 
the set of conceivable future events of which at least one must be realized, and proceeding 
to the one that is first to happen, i.e. the one with the largest t among them. That is why 
we must be able to calculate the value of t for each of those events. 

For all events we will give a polynomial equation describing it. The equations will be 
solved either for t or Oi(t). Once we have Oi(t) we can easily obtain t from the following 
inverse function of Oi(t). 

2( 0i (t) -mi) v ' 

In the process of determining t, we will repeatedly encounter fixed degree polynomial 
equations. Solving them can be assumed to be a constant time operation, see [2]. 



4.2.1 Jumping and retargeting 

Jumping is any event in which C(t) abruptly changes its position by switching to some 
other hyperbola segment. In order to efficiently find state transition events that include 
jumps, we always need to know which hyperbola segment can be jumped to from the 
current position. There is always at most one such target hyperbola segment, and we will 
see how to keep track of it. 

Consider some point H{x) on the plot. Let jump(x), the jump destination for x, be the 
index of a hyperbola segment which contains the rightmost point H(w) on the plot such 
that w < x and the ray from H{x) through H(w) only touches the plot, but does not 
intersect it properly. That is, jump(x) is the index of the lowest visible hyperbola segment 
when looking from the point H{x) to the left. If there is no such w, because hyperbola 
segments on the left are obscured by the one containing H(x), then we set jump(x) to be 
the index of the hyperbola segment containing H(x). 

Consider only the values of x where jump(x) changes. We call such values retargeting 
positions, and points H(x) retargeting points (Figure Yty. 



Figure 7: Retargeting point. 

Theorem 1. It is possible to find all retargeting positions, ordered from left to right, 
together with jump destinations of those positions in 0(n) time. 

Proof. The algorithm for finding all retargeting points is similar to finding the lower convex 
chain in Andrew's monotone chain convex hull algorithm [3], except the sorting part is 
not needed. 
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Let jo be the index of the hyperbola segment which contains any of the lowest points of 
the plot of h. We scan segments from left to right and maintain the longest chain of those 
hyperbola segments that are in a convex position, starting with the T~Lj and ending with 
the segment currently being processed. Intuitively, convex position means that for each 
three hyperbola segments appearing consequently in the convex chain, left and right of 
those segments cannot be connected with a line segment passing completely under the 
middle one, see Figure [8j 



Figure 8: Hyperbola segments drawn with a solid line are in convex position. Remember 
that each hyperbola segment contains its left endpoint, but not the right one. 

The stack is used for keeping the track of segments in the chain (we store their indices 
only). At each of the following segments encountered, we pop segments from the stack 
until the current segment can be added to the chain so that the chain is again in convex 
position, after which we push the current segment onto stack. 

Algorithm 1 Retargeting Points 
RetargetingPoints [ 
Push(jo - 1) 
Push(jo) 

for k <— jo + 1 to jo + n do 
loop 

i <— Stack second 
j <— Stack top 
li «— Tangent (i, j) 
h «— Tangent(j, k) 
if l x n U k 7^ then 

Let g be the leftmost point from l\ n %k- 



if l 2 C\U k / then 

Let g be the leftmost point from I2 n %k- 

Append g to RetargetingPoints, and set jump destination of g to j. 
Push(/c) 
end if 
break loop 
end if 
end loop 
end for 



Tangent (i, j) returns the line which acts as a boundary between the region of the plane 
from which the Hi is visible from below of T~Lj and the region from which Hi is obscured 
by the T~Lj. More precisely, TANGENT(i, j) is the line touching both cl("Hj) and cl('Hj) from 
below, where cl(X) denotes the closure of a point set X, which is, in our case, essentially 
an addition of the right endpoint to the hyperbola segment. 





else 



S 



Each point reported is indeed a retargeting point since the jump destination changes at it. 
Points reported in the "if branch have the property that points on Hk immediately left 
and right from g have Hj and Hi as their jump destinations, respectively. Points reported 
in the "else" branch have the property that points on Hk immediately left and right from 
g have Hj and Hk as their jump destinations, respectively. 

Let us make sure that no retargeting points were omitted. Each retargeting point g on 
Hk must lie either on Tangent^, j) or Tangent(j, k), for some i < j < k. In the first 
case, Hi and Hj must be two consecutive elements of the convex chain ending with Hk-i, 
otherwise Hi or Hj would not be visible from the neighbourhood of g. Since Tangent^, j) 
intersects Hk, the same must be true for each subsequent pair from the convex chain on 
the right of i and j. Thus, while the algorithm pops elements from the stack, in one 
moment the last two elements on the stack will happen to be exactly i and j. At that 
iteration of the inner loop g will be reported in the "if branch as the retargeting point. 
The other case, when g is on the Tangent(j, k) is similar: Hj must be the element of the 
convex chain ending with Hk-i, otherwise it would not be visible from the neighbourhood 
of g. The tangent of each subsequent pair from the convex chain with the right element 
on the right of j intersects Hk- Thus, the algorithm will be popping elements from the 
stack until the top element of the stack is j, when g will be reported in the "else" branch 
as the retargeting point. 

The running time of algorithm is O(n), since TangentQ can be computed in constant 
time, and each hyperbola segment index k £ {jo + 1, . . . , jo + n} is pushed on stack and 
popped from stack at most once. Hence, the number of retargeting points reported is also 
O(n). 

Retargeting points reported by the algorithm come in order sorted from left to right. This 
is because the retargeting points reported in a single iteration of the outer for- loop belong 
to the same hyperbola segment, and segments come in left-to-right order. Retargeting 
points reported on the same hyperbola segment are also in the left-to-right order: inside 
the inner loop, Hk is consecutively intersected with lines such that each line is of the lower 
slope than previous and lies beneath it under Hk- Hence, each subsequent intersection 
point lies to the right of the previous one. 

Here we reported only retargeting points from a single period of the plotted function, but 
all other can be found by simply translating these horizontally by the integral number of 
periods fi(P). □ 

Next, we list all types of events that can happen while moving the lever leftwards and 
show how to calculate the value of t for each. When talking about events that involve 
jumps, the current jump destination, jump(c(t)), will be denoted by j m . 

4.2.2 Jump destination change event 

Jump destination j m changes whenever C(t) passes over some retargeting point. At that 
moment it is necessary to recalculate all future events which involve jumps. Let z be the 
next retargeting position, i.e. the rightmost one that lie to the left of t c . Depending on 
the lever state, we calculate the event position in one of the following ways. 
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• If the current state is — The jump destination change event cannot occur 
before leaving this state since C(t) stands still at the "wedge tip", so it cannot pass 
any other point. 

• If the current state is ~ The lever is tangent to the ~Hj, so this event can only 
happen if z > rrij. Otherwise, the lever would have nonpositive slope when touching 
the plot at H(z). The equation describing this event is 

Cj(Oi(t)) = z, 

and it solves to 

d 2 — zm-j + m 2 

Oi(t) = - 1 

raj — z 

• If the current state is - Right endpoint of the lever slides over %j and will 
coincide with H(z) when 

t = z-n(P)/2. 

Note that j m is not used to describe lever state, so jump destination change event does 
not change the current state. All other considered events, however, are state transition 
events. 



4.2.3 ^ 

The interval to which t belongs changes from [e^, ej+i) to [ej_i, ei). This event happens at 
ej. 

t = ei. 



4.2.4 (i,j)y -> (ij-l) y 

This is the event when the right endpoint of the lever simply slides from one hyperbola 
segment to another. 

t = ej- fi(P)/2. 

4.2.5 {i,j)y -> (i, j) K and (i, j) K -> (ijf 

In this event the lever changes from being a tangent to Hj to touching Hj with its right 
endpoint, or the other way round. The corresponding equation for this event is 

c j (o i (t))=t + ti (P)/2, 

which can be transformed to a cubic equation in t. Since there can be at most three real 
solutions to that equation, it is possible that this event takes place at most three times 
with the same % and j. On each occurrence of the event the lever switches between being 
a tangent and touching the plot with its right endpoint. 



10 



4.2.6 (i,j) y ^(i,j m ) K 

This event happens when the lever state changes from having an endpoint on Hj to being 
a tangent to %j m . The equation is 

* 3m W))- h . {t)+fl{ py 2 > 
which further transforms into a 5th order polynomial equation in t. 

The line through Oj(i) with slope Sj m (oi(t)) touches the hyperbola hj m , but we need to be 
sure that it actually touches the segment %j rn of that hyperbola. It may as well be the 
case that Hj m is not wide enough to have a common point with the line. More precisely, 
the first coordinate, u, of the touching point between the line and hj m must belong to 
interval [e Jm , e Jm+ i). To get that coordinate, we solve the equation 

h j m ( u ) = S jm(°i(t))- 

Considering that Oi(t) < rrij m < u must hold, we get a single solution 

u = m jm + 



m jrn -Oi(t)' 



If u ^ [e Jm , e Jm +i), we do not consider this event. 

This check, to see if the line through Oi(t) with slope Sj m (oj(t)) actually touches the 
hyperbola segment Hj m , we call collision check, and we will also use it in some other 
events. 



4.2.7 (i,j) y ^{i,j m ) V 

The event when the lever state changes from having an endpoint on T-Lj to touching the 
wedge between %j m _ 1 and Hj m is described by 

h jm {e jm ) fy(t + /i(P)/2) 



e jm -Oi(t) hi{t)+^P)/r 
which again transforms into a 5th order polynomial equation in t. 

4.2.8 (i,j) K ^(i,j) v 

This event happens when the point in which the lever is touching T-Lj reaches ej. Here, 
the lever is tangent to Hj, and since it must have a positive slope, this will only happen 
if ej > nij. The event equation is 

Cj( 0i (t)) = ej, 

which solves to 

d 2 , - ejrrij + m\ 

o l {t) = -± — >-. 

m, — e, 
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4.2.9 (i,j) K -> (i,jm) K 

This event happens when the lever becomes a tangent to two hyperbola segments, T~Lj and 
%j m simultaneously. It can only happen if %j m is lower than Tij, i.e. dj m < dj. 

S jm(Oi(t)) = Sj(Oi(t)). 

Since Oi{t) < rrij m and Oj(i) < rrij, the only solution is 

lj ~ d--d- 



Here we need to apply the collision check explained in Section 4.2.6 to see if the common 
tangent actually touches %j m - If the test fails, we do not consider this event. 



4.2.10 (i,jf -> (i,j m ) v 



The event in which the lever touches the wedge tip at ej m while being a tangent to hj is 
represented by the following equation. This can only happen if hj m (ej m ) < dj. 

himieSm L = 

Cj m -Oi(t) 

This can be transformed to a quadratic equation in Oj(i). The two solutions correspond to 
two tangents to hj from the point (e,- m , h(ej m )). The smaller of the two solutions is where 
this event happens. 



4.2.11 {i,j) v ^ (ij-l) y 

This event happens when the lever stops touching the tip of the wedge and starts to slide 
its right endpoint over the hyperbola segment on the left of the wedge. 



t = e j - t x(P)/2. 



4.2.12 



X 



This event happens when the lever stops touching the tip of the wedge and becomes a 
tangent of the hyperbola segment on the left of the wedge. This can only happen if 
ej > rrij-i. 

Cj-i(oi(t)) = ej, 

which solves to 

mj-i- ej 
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4.2.13 (i,j) v ^ (i,j m ) 



This event happens when the lever stops touching the tip of the wedge and becomes a 
tangent of the hyperbola segment %j m . This can only happen if hj(ej) > dj m . 



From that we get a quadratic equation in Oi(t). The two solutions correspond to two 
tangents to hj m from the point (ej,h(ej)). The smaller of the two solutions is where this 
event happens. 



4.2.14 (i,j) v ^ (i,j m ) V 

This event happens when the lever touches two wedges, at points Ej m and Ej simultane- 
ously. The condition for that is 

ej-Oiit) e jm -Oi(ty 
which is a simple linear equation in Oj(i). 



4.2.15 Sequence of states 

We want to find a sequence of states through which the lever will pass on its leftward jour- 
ney, together with the positions where state changes happen. Let the obtained sequence 
be pi,Si,p2,S2,P3, ■ ■ ■ ,p r ,<S r , where p\ < p2 < . • • < p r - Each state takes part when 
the lever position is exactly between pk and Pk+i, where p r +i = Pi + m(-P). We call this 
the sequence of realized states. 

To calculate the sequence of realized states, we start from a specific lever position that 
has a known state. Let pi ow be any of the values for which h(pi ow ) is minimal, and 
T~Lj be the hyperbola segment containing it. The algorithm starts with a lever in position 
t c = to = Plow — A t (-P)/2- This lever has its right endpoint on the plot at the point H{pi ow ), 
which means that its state is (io,jo)^, where iq is the index of hyperbola segment over to. 
pi ow is also a retargeting position, so we also know initial jump destination. 

The algorithm then iterates with the following operations in its main loop. It first calcu- 
lates all possible events that could happen while in the current state. Among those events 
let E be the one with the largest t that is not larger than current t c . It is the event that 
must occur next. The algorithm sets t c to equal t, and updates jump destination if E is 
jump destination change event, or switches to the new state if the event is state transition 
event. In the latter case, position t and the new state are added to the sequence of realized 
states. These operations are iterated until one whole period of the plot is swept, ending 
with t c = to — fi(P)/2 in state (io — n,jo — n)^ . 

Theorem 2. Sliding Lever Algorithm runs in 0(n) time, and the length of the 
produced sequence of realized states is 0{n). 
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Algorithm 2 Sliding Lever Algorithm 



Find pi ow and j . 

to <- Plow ~ M( p )/ 2 
Find io- 

Run Retargeting Points to find retargeting points and their jump destinations. 

t c <— to 
i 4- i 

j <- jo 

jm ^— jump destination of pi ow 
Set the current state to (io,jo) ■ 
Add to an d (iojio) to the sequence, 
while t c > to — n{P) do 

Calculate all the events for the current state. Ignore jumping events if j = j m . 

Let E be the first event to happen (the one with the largest t not larger than t c ). 

t c <— t of the event E. 

if E is jump destination change event then 

Update jm- 
else 

Set the current state to the destination state of E. 
Add t and the current state to the sequence, 
if E is a jumping event then 

j m <~ j 

end if 
end if 
end while 



Proof. Each event is either jump destination change event or state transition event. From 
Theorem [T] we have that there are 0(n) jump destination change events, and now we will 
show that there are O(n) state transition events. 

Each state transition event transitioning from some {i,j) state decreases either i, or j or 
both. The only exception are the events (i, j) —> (i,j)^ and (i,j) — > (i,j) , however 
those events can happen at most three times in total for the same i and j. Note that 
jm < j, but when j m = j, we do not consider events involving j m . Variables i and j start 
with values io and jo, and, after the loop finishes, they are decreased to iq — n and jo — n. 
Hence, no more than 0{n) state transition events occurred, implying the linear length of 
the sequence of realized states. 

Calculation of each state transition event takes a constant time, there is a constant number 
of events considered, and loop is iterated 0(n) times. To find next jump destination change 
event, we move through the sorted list of retargeting positions until we find the first 
retargeting position not greater than t c . The total time for calculating jump destination 
change events, over all iterations, is linear. Finally, the running time of the whole algorithm 
is also linear. □ 

While choosing the next event, we did not consider the possibility that there are several 
events with the same largest t not larger than t c - But, if that happens, we choose an 
arbitrary one to be the next event. All of them will be realized in some form while lever 
position equals t, and no matter the order we processed them, the lever will end up in the 
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same state once t c passes t. 



5 Merging the two dilations 



Knowing the sequence of realized states is sufficient to determine the exact lever slope at 
any position. Remember, the lever slope at position t is the inverse of the left dilation via 
P(t) Q. But, to know the dilation via some point we need both left and right dilations 
via that point ([!]) . 

Sliding lever algorithm was designed only for left dilation, but an analogous algorithm 
can be designed in the same manner for right dilation (or by performing the algorithm 
for left dilation on the mirror image of the polygon P, and then transforming obtained 
results appropriately). This implies the concept of the right dilation lever for t (as opposed 
to the left dilation lever, or just lever, as we have been calling it until now), which has 
negative slope and touches the plot on the left side of t. This is where we return to the 
original notation of + and — in superscript denoting relation with left and right dilation, 
respectively. 

Let pt,S^,pt,S^,p^, . . . and pJ,S{ ,p% ,p% , ■ ■ ■ ,P~-,S~_ be the sequences 

of realized states for left and right dilation, respectively, where both sequences p + and p~ 
are in nondecreasing order. States for right dilation are described by notation as 

well, with the meaning analogous to the meaning of the notation for left dilation states. 
We say that the (right dilation) lever for t is in the state when Hi is the hyperbola 

segment above t, and the (right dilation) lever touches the hyperbola segment tij- 

Having obtained both sequences of realized states, we can now merge them into a single 
sequence pi,«Si,p2>«->2)P3> • • • ,Pr,Sn where p\ < P2 < . . . < p r is sorted union of p + and 
p~ sequences. States of the merged sequence are pairs consisting of S + and S~ state. Let 
p^ + be the last position from the left dilation sequence appearing before state S in the 
merged sequence, and correspondingly, let p7 be the last position from the right dilation 
sequence appearing before state S in the merged sequence. State S is then ordered pair 

By the Theorem [2J both r + and r~ are 0(n), so the length of the merged sequenced is 
also linear in n. Hence, the merged sequence can be computed by a simple 0(n) time 
algorithm. 

For each state Sk = > S7_ ) we have a single expression for computing the lever slope 
when pk < t < Pk+i, both for the left dilation and for the right dilation. To find minimal 
dilation while in that state, we want to find t such that the minimum of the two slopes 
for the two levers for t (left and right) is maximal (fil pi) 



mm op( t \ = ; — -r-T-, 101 

Pk<t<Pk+i max Pk<t<Pk+i mm{s+(t), s (t)\ 

where s + (t) is the same as s(t), the slope for the left dilation lever for t, and s~(t) is the 
slope for the right dilation lever for t. 

Let us analyse the shape of the functions s + (t) and s~ (t). Assume that the corresponding 
state where t belongs is S = (S + ,S~). 
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If S + is a (ijj) 1 ^ state, then, from equation (|4j), we have 



hj(cj(oi(t))) f m j ~ Oi(t) 
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Lemma 3. If S + is a {i,]) 1 ^ state, then s + (t) is monotonously increasing function. 

Proof. From ([7]) we see that Sj(oi(t)) is monotonously increasing when Oi(t) < mj which 
holds in the given state. Since o(t) is monotonously nondecreasing (Lemma [2]), it means 
that their combination, s + (t) is monotonously increasing as well. □ 

If S + is a state, then we have 

8+{t) = 8+ (t) = J^L. (8) 

V>3) ej-Oi(t) 

Lemma 4. If S + is a (i,j) v state, then s + (t) is monotonously increasing function. 



Proof. From ([8j) we see that s + (t) is monotonously increasing in terms of Oj(i) when 
Oi{t) < ej which holds in the given state. Since o(t) is monotonously nondecreasing 
(Lemma 0), it means that s + (t) is monotonously increasing in terms of t as well. □ 

If S + is a (i,j)^ state, then we have 

* (*)-^p(*)-^. (t) + M(p)/2 > (9) 
which is not necessarily monotonous. 

Similar observations hold for right dilation analogues: s~ (t) is monotonously decreasing 
if S~(t) is (i,j) lc or (i,j) v state. 

Lemma 5. is a (i,j) y state then s + (t) < s~(t). 

Proof. From equation ([2]), using the fact that hj(t) = hj+ n (t + fj,(P)) holds because of the 
periodicity of the plot, we have 

> min s(t,t r ) = s + (t). 
t r e[t,t+n(P)/2] 

□ 

Analogously, if S + is a (i, j) y state then s~(t) < s + (t). 

We need to calculate max Pfc < 4 < Pfc+1 min{s + (i), needed in equation (|6j), which is 

equivalent of finding the highest point of the lower envelope of the functions s + (t) and 
s~(t) (Figure^). This calculation depends on what the types of states S + and S~ are. 
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Figure 9: m&x Pk < t < Pk+1 mm{s + (t) , s (t)} 



S + is (-, state, 5 is (-, state : Using LemmaJHJwe get s + (t) = s (i), so 

max min{s + (t), s _ (i)) = max s~!~. ., v (t). 

Pk <t< Pk+1 1 W ' W/ Pfc <i<P fe+ i <«>* V ' 

The maximum is achieved either at interval ends or at local maxima, which are 
obtained by solving a polynomial equation. 

S + is (•, state, S~ is (•, state : Using Lemma[5]and Lemma[3]we get 

max min{s + (t),s"(t)} = max sj )K (t) = s+ j)K (pk+i)- 
Pk<t<Pk+i Pk<t<Pk+i W 

5 + is (•,-)' c state, 5~ is (•,-)' c state : From Lemma |3j we know that s+_. + ^(i) is 
monotonously increasing, and sT. is monotonously decreasing. The highest 

point of the lower envelope of their plot on \pk,Pk+i] is thus located either at one of 
the interval endpoints, or at the point of the intersection of two plots, which can be 
found by solving a polynomial equation. 



Other six combinations of state types, which are not listed, are resolved in the same 
manner as the three mentioned above. 

Finally, by taking the minimal of all r dilation minimums from \pk,pk+i\ intervals we 
obtain the overall minimum dilation. 

5 = min min Sp( t \. 

ke{l,2,...,r} P k<t<Pk+i w 



While calculating the minimums we maintain the value of t for which the minimum is 
achieved, so we also get the point on P which is the endpoint of the optimal feed-link. 
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